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GENERAL ALGEBRAIC METHOD APPLIED TO CONTROL ANALYSIS OF 

COMPLEX ENGINE TYPES 

By Aaron' S. Boksenboh and Richard Hood 


SUMMARY 

A general algebraic method of attack on the problem of con- 
trolling gas-turbine engines having any number of independent 
variables was utilized employing operational functions to 
describe the assumed linear characteristics for the engine, the 
control, and the other units in the system. Matnces were used 
to describe the various units of the system, to form a combined 
system showing all effects, and to form a single condensed 
matrix showing the principal effects. This method directly led 
to the conditions on the control system for noninteraction so that 
any setting disturbance would affect only its corresponding con- 
trolled variable. The response-action characteristics were 
expressed in terms of the control system and the engine char- 
acteristics. The ideal control-system characteristics were ex- 
plicitly determined in terms of any desired response action. 

INTRODUCTION 

The current- development of gas-turbine engines indicates a 
future trend toward a wide variety of engine types. As new 
engines are developed or built by combining basic components 
of existing engines, the control problem presented by each 
engine type will be different. Control systems satisfactoiy 
for one engine type could not be expected to operate for 
another engine type. Even engines of the same type but of 
different size or operated under radically different conditions 
may have different engine characteristics that require special 
control characteristics. The current fluid and indeterminate 
state of gas-turbine development indicates a need for a gen- 
eral control analysis that could be applied to any of the new 
engine types as they are developed. 

The control problem for engines' differs from that for many 
plants, servomechanisms, or electric networks, because the 
control of more than one variable is desired. At least those 
variables that may cause failure should be controlled. For 
gas-turbine engines, failure may be caused by excessive tem- 
perature, speed, or torque. Exceeding specified burner oper- 
ating limits may cause burner blow-out. The problem of 
controlling more than one variable is therefore basic for gas 
turbines. 

In control design, limitations may be placed on such char- 
acteristics as response time, maximum undershoot or over- 
shoot, dead-band, or steady-state error. This problem is 
the subject of much of the control literature to date and 
various methods of exact, approximate, and graphical means 
of solution have been developed (reference 1). 


Controlling more than one variable generally introduces an 
interaction among the controlled variables. Making a new 
setting for one variable may cause, during the transient 
state, changes in the other controlled variables. If these 
other variables are operating at or near a maximum point, 
this interaction may cause excessive values and possible 
damage to the engine. It would therefore be desirable that 
each new setting of a controlled variable affect only the 
variable being set, thus giving separate noninteracting con- 
trol of all the variables being controlled. 

A general algebraic approach to control problems of multi- 
variable engines, which was developed at the NACA Lewis 
laboratory during 1949, is presented herein. This method 
leads to the conditions on the control characteristics that 
eliminate interaction between control variables and to an 
explicit solution for control characteristics that yields any 
desired response action. 


SYMBOLS 


The following symbols are used in this report: 
a, b, c, e constants 

C or C (p) control function to which engine-dependent- 
variable. errors are applied 

C' or C' ip) control function to which engine-independent- 
variable errors are applied 
E or E ( p ) engine-characteristic function 
E* square matrix using first i columns of engine 

matrix 

\E*\ determinant of E* matrix 

\E* Jt \ cofactor of E Jt element of E* matrix 

Ftp) general operational function 

i number of engine-dependent variables being 

controlled 

L or L ip) instrument-characteristic function 
N engine speed 

N . measured value of engine speed 

n number of engine-independent variables 

p operational symbol 

R or R ip) response function for controlled engine- 

dependent variables 

R' or R' (p) response function for controlled engine- 

independent variables 
S or S (p) servo-characteristic function 
T engine temperature 

T measured value of engine temperature 
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uncontrolled transient disturbance in engine- 
independent variable 
engine fuel flow 

measured value of engine fuel flow 

signal to engine fuel-flow servo 

setting of engine-independent variable 

engine-independent variable 

measured value of engine-independent variable 

signal to engine-independent variable servo 

setting of engine-dependent variable 

engine-dependent variable 

measured value of engine-dependent variable 

propeller-blade angle 

signal to propeller-blade-angle servo 


engine-interaction factor 



engine time constant 


indices 

setting of specific engine variable 


In addition to the intentional disturbance made on the 
control setting, inadvertent transient disturbances may 
appear in any part of the system. One possible source of 
such transient disturbances has been included as an uncon- 
trolled change in engine-independent variable. The variable 
x is thus the sum of a controlled part and of an uncontrolled 
part F. 

The characteristic functions used to describe the various 
units of the system are operational functions or Laplace 
transforms. The dynamic characteristics of gas-turbine 
engineSj _ controls, and servos can be approximated by such 
operational functions at least for the small deviations of 
variables in which the linear assumption is valid. These 
operational functions can be used without reference to the 
initial conditions of operation before disturbances occur, if 
the initial condition is a steady-state condition and the values 
of all variables are their deviations from the initial steady- 
state condition. In general, the unit can be described by a 
general operational or transfer function as follows: 


y—F(p) ■ x 


( 1 ) 



ANALYSIS 

The system considered herein is essentially that shown in 
figure 1. This block diagram describes a typical engine and 
control system. The general system configuration chosen 
is the one employing negative feedback in which error in the 
controlled variable is fed back into the control. The engine 
is represented as a box with a characteristic function E, the 
input of which is the engine-independent variable x and the 
output is the engine-dependent variable y. The characteris- 
tics of the instrument measuring the value of the controlled 
variable are denoted as L for which the input is the variable 
itself and the output is the measured value y. The error is 
the difference between the setting Y and the measured value 
of the variable y. The characteristics of the control are 
denoted as a function C for which the input is the error and 
the output is a signal x to a servo unit, which is the device 
that actually changes the engine-independent variable and 
completes the loop. The servo characteristic is. denoted as 
the function <S and may be considered as a part of the control 
or as a separate entity, as shown in figure 1 . 


The block diagram for this unit can be shown as a box 
containing the function F(p ) for which the input is x and the 
output is y. 

One method of deriving the operational function for a unit 
that approximately follows a linear differential equation is as 
follows for the equation given by 

y+ay+by=cx+ex (2) 

where the dots above symbols indicate derivatives with 
respect to time. This equation can be algebraically repre- 
sented in terms of the differential operator or in terms of the 
Laplace transform variable (reference 2). If the differential 
operator p is substituted for d/dt, equation (2) can be written 

p 2 y+apy+by=cpx+ex (3) 

If this equation is solved for y, an explicit expression for the 
output can then be obtained as a function of p operating on 
the input x as follows: 

(4) 

For each unit of the system, an equation similar to equa- 
tion (4) can be written in terms of the input, the transfer func- 
tion, and the output. A set of simultaneous algebraic equa- 
tions is then obtained that can be combined or manipulated in 
any desired manner. The resulting transfer functions can then 
be interpreted from the domain of the variable p to the real 
time domain as time responses of a variable to some specified 
disturbance. 
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For the system of figure 1, the following equations can be 
written: 

Engine characteristic: 

y=E- x 

Control characteristic: 

x=C-(Y-y) 

r (5) 

Instrument characteristic: 

y=L-y 

Servo and transient disturbance: 
r=fS-x+I ;r 


The functions E }t represent engine characteristics whereby 
each input separately affects each output. 

Equation (7) or equation (8) can be visualized in the matrix 
form of figure 2. The engine functions are set up in a rec- 
tangular array where E jt represents the element in the jth 
row and the kih column. The x inputs enter the matrix in 
the columns and the y outputs are attached to the rows. 
Each input multiplies all the elements in its column and each 
output is the linear sum of the resulting products in its row. 
Thus the effects shown by equations (7) and (8) can be 
pictorially traced in figure 2. For later use, a square matrix 
obtained by using only the first i columns in figure 2 is 
denoted by E*. 


Equations (5) can be combined for the effect of the disturb- 
ances Y and V on the controlled variable y as follows: 


ESC , E 
y ~ESCL+l' 1 ' ESCL+l' 


( 6 ) 


The previous procedure is, in general, followed herein. 
Equations, such as equations (5), are written for each unit 
in the system and are then combined, as in equation (6), to 
obtain the effects of the disturbances. 

General representation of engine. — If the engine has more 
than one independent variable, there are additional degrees 
of freedom in control. The engine-characteristic function 
must be expanded to show how each dependent variable of 
interest is affected by all the engine inputs as follows: 


yi=Eu • Xi-\~Ei3 • x 3 + . . . -\-Ei„ • x,' 

y 2 =En • x t +J ?22 • x 2 + . . . +f? 2 * • x„ 

. . 


(7) 


y { =En • Xi+E< 2 ■ z 3 + . • ■ -\-E tx • x„ J 


In equations (7) each dependent variable y is expressed as a 
linear sum of effects due to each independent variable x. For 
n independent variables, there are n terms in each equation. 
There is an equation for each dependent variable of interest. 
For i dependent variables, where i^n, there are i equations. 

The number of degrees of freedom of the entire system 
cannot be greater than the number of degrees of freedom of 
the engine alone. For continuously acting controls (excluding 
such controls as limiters) , the total number of engine variables 
being controlled should not be greater than the number of 
engine-independent variables. Equations (7) can be written 
in the following index form: 


y s — Eft • Xt ( 8 ) 

t=i 

where 

j == 1,2, . . . % 


Xt x 3 x n 



General representation of control. — The control system is 
generalized to this extent: (1) Inasmuch as the engine has n 
independent variables, there are n degrees of freedom for the 
entire system. If only a limited number i of engine-dependent 
variables are to be controlled, the difference n—i of the 
engine-independent variables can be controlled. (2) The 
control system employs negative feedback where only errors 
are applied to the control. (3) Each error is to affect 
every engine-independent variable. 

The control equations are written where each control 
output x is expressed as a linear sum of effects due to the 
errors in the controlled variables as follows: 


x l = C n -(Y l -yi)+C 12 -(Y i -yJ+ . . . C u (Y t -y t )+ 

• C' Ki+u • C-Xi+i r 4+ i)+ ~\~C i*(A„ x») 

x*=C n i-(Yi— yi) + C«2-(F 2 — yi)-\- - - - +C n r{Y i —yi) J r 

1) ’ an. ’ C*!* x n) * 
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There are a total of n controlled quantities, n errors, and n 
terms in each equation. There are also n engine-independent 
variables and therefore n equations. Those independent 
variables ta.be controlled are given subscript numbers from 
i-fl to n and those controls to which the corresponding 
errors are applied are denoted by C . . 

The control characteristics of equations (9) can be written 
■ in the following equation in index form; 


output, the measured value. These characteristics can be 
written as 

y,=L e , ■ y, 

where . 

o = l,2, . . . i 


where 


Xjx — ■ Xu 

n=i+l, . . . nj 


( 11 ) 


&=SCi. * (Y„—y l )+ S • (.Xp—Xp) ' (10) 

«=i p-i+i 

where 

k= 1, 2, ... n 

Equations (9) and (10) can be visualized in the matrix form 
of figure 3. The control functions are set up in a square 
arrangement where C k , represents the element in the &th 
rov r and oth column. The error inputs, Y—y and X— x, 
enter the matrix in the columns and the x outputs are 


The general assumption is made that each engine- 
independent variable would require a servo, which would 
change the variable in accordance, with its corresponding sig- 
nal x. A source of transient disturbance is included that may 
vary any of the engine-independent variables outside of the 
controlled variation in them. Thus, the expressions for the 
engine-independent variables can be given by 

Xt=Sici • £*+ T r * (12) 



attached to the rows. Each input multiplies all the ele- 
ments in its column and each output is the linear sum of the 
resulting products in its rov r . Thus, the effects shown by 
equations (9) and (10) can be pictorially traced in figure 3. 

Representation for instrument and servo functions. — The 
characteristics of the instrument measuring the value of the 
controlled variable have been denoted as L in equations (5) ; 
the input of this instrument is the variable itself and the 


Each variable x is thus the sum of a controlled part and an 
uncontrolled part. 

Analysis of complete system. — The determining charac- 
teristic equations for all the units in the system are repeated. 


Engine: 


where 


y } =±E ik -x k 
.7 = 1,2, . . . i 


( 8 ) 


Control: 


& CV( Xp-%) (10) 

5=1 P-i+1 

where 

k = 1,2, ... n 


where 


where 


where 


Instrument: 

o=l,2, . . . i 

Xp — * Xp 

n=i+ 1, . . . nj 
Servo and transient: 

Zk= S M ■ x k + V k 
k= 1,2, . . . n 


^ (13) 


(H) 


( 12 ) 


These equations can be visualized in the matrix form of 
figure 4. Figure 4 represents, by block diagram, the general 
system considered herein. This figure specifically is for an 
engine having five independent variables for which three 
engine-dependent and two engine-independent variables arc 
controlled. The engine matrix is shown for the five inputs 
and the three outputs to be controlled. The control matrix 
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Transients 



Figure 4. — Combined general eontxoUed -engine-system configuration. 


is shown for the three y error inputs, the two x error inputs, 
and the five outputs to the engine servos. The measuring 
devices, the servos, and the transients are also included. 
The entire system is enclosed except for the outside dis- 
turbances, which can be imposed on the system. These are 
the five settings of controlled variables T and X and the 
five transient disturbances that have been included. 

All the properties of the system can be obtained from this 
block-diagram representation. The effect of any outside 
disturbance or any internal disturbance can be easily traced 
on figure 4. For specific effects, the equations of each unit 
of the system, equations (13), are most productive of results. 

When equations (13) are combined to give the effects on 
the controlled variables of the outside disturbances, the 
following equations are obtained: 

2/j=2j S Ejk&uCttil , — L„y,)-{- 

*=I V = 1 

S X EgtVt (14) 

t = l jt = i+l Jt =1 

where 

3= 1.2, . . . i 

x k = S S a C k ,( Y,- L„y ,) + ± SuC'^X,- L m x,) + V k (1 5) 

v-l (i=i+l 


Equation (15) is valid for the effects on any (k=l, 2, . . . n) 
engine-independent variable. For the effects on the con- 
trolled engine-independent variables, equation (15) is used for 
fc=i-f 1 to ii. A more condensed form of the matrix repre- 
sentation of the system, shown in figure 5, is indicated by' 


r, x t 



A J 

r 

A 

12 


n, T 1 

1^ 
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Vt 


Vi 
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System matrix 


Fig CEE 5.— Condensed matrix form of complete controlled system. 


equations (14) and (15). A single-system matrix is used in 
which the inputs are the errors in controlled variables and 
the outputs are the controlled variables. The ESC matrix 
is a matrix in which the element in the jth row and rth 

n 

column is IY.E^Sh-Ch,. Similarly, an ESC ' matrix is a 

t=i 

matrix In which the element in the yth row and pith column 

n 

is y'.EttSttC' to. An SC matrix is defined as S^Ct, and an 

*=i 

SC' matrix is defined as SaC'^. An additional matrix is 
included showing the effects of the transient disturbances. 

Noninteraction conditions. — Several types of noninter- 
action may be specified as follows: 

(1) A specific controlled variable, such as y } , is to be 
affected only by its own setting and is not to be affected 
by any other setting, which means that in figure 5 the jth 
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row of the system matrix contains only the diagonal element 
and the other elements must be zero. 

(2) Any specific setting, such as Y tl is to affect only its 
corresponding controlled variable y, and not any of the other 
controlled variables. In figure 5, the Zth column of the 
system matrix contains only the diagonal element and the 
other elements must be zero. 

(3) Every controlled variable is to be affected by its 

corresponding setting only, or every setting is to affect its 
corresponding controlled variable only. This complete non- 
interaction means that in figure 5 the system matrix is a 
diagonal matrix. In general, the effect of the setting on the 
system is of primary interest and a condition such as (1) is 
of secondary interest. The interaction effects defined in 
terms of the effects of the settings, such as (2) and (3), will 
thus be considered herein. . 

In order for any setting such as Y ( to affect only y t , it is 
sufficient from figure 5 that for all j from 1 to i when jv&t 


for all j from 1 to i; and for all y from t+1 to n where y^r 

S llfL C'„=0 (21) 

or 

C '„= 0 (21a) 

Equation (21) states that the i+1 to n elements of the rth 
column of the O' matrix contain only the diagonal term. 
With this condition, equation (20) becomes 

(.Ej k S kk O > k ^-^- Ej r S rr C' r r=0 (22) 

for ally from 1 to i. Equation (22) expresses i homogeneous 

equations in the i+1 unknown (C')’s and therefore specifics 
the ratios between any two elements in the rth column of 
the C matrix. This result becomes 


X Ej k S kk Cm — 0 (16) 

*=. x 

and for all y from i+1 to n 

5^=0 (17) 

or 

<?„ ( =0 ' (17a) 

Equation (17) states that the i+1 to n elements of the ith 
column of the C matrix (fig. 3) are to be zeroed. Along with 
this condition, equation (16) becomes 

X! E Jk S kk C kt =0 (18) 

for all j from 1 to i where j^i. Equation (18) expresses 
i— 1 homogeneous equations in the i unknown C’s and there- 
fore specifies the ratios between any two elements in the ith 
column of the C matrix. This result becomes 


EjjCj ( _ 

E\j 

&vv(svt 

E\, 


(19) 


as shown in the appendix. The notation E* refers to the 
square engine matrix using only the first i columns. The 
notation \E* tj \ indicates the cofactor of the E tj element’of 
the E* matrix. The necessary and sufficient conditions for 
any setting Y t to affect only y t are 

C„ t = 0 (17 a) 

for all y from i+1 to n and 


as shown in the appendix. 

For any setting X, to affect only x r , it is sufficient from 
figure 5 that 


&! 

E*J 


(19) 


q pi Sl-E*wl^*r 
OjjO Jr k-l 

nr — i 


(23) 


for any j and v from 1 to i or 


EpC’jr 

SrrC'„ 



(24) 


as shown in the appendix. The necessary and sufficient 
condition for any setting X r to affect only x, is 


CV= 0 

for all y from i+1 to n when y^r and 
or Q / — '52l\E*tj\E kr 

Ojjly Jr It - 1 

SrrC'rr |S*I 


(21a) 


(24) 


as shown in the appendix. 

For complete noninteraction, where every setting i3 to af- 
fect its corresponding variable only, the conditions can be 
summarized as follows: 


<?*= 0 (17a) 

where 

y=i+ 1, . . . n 
t= 1,2, . . . i 

<?V=0 (21a) 

where 

M=i+ 1, ...» (y^r) 

r=i+l , . . . n 


XJ Ej k SnC'to — 0 

fc= 1 


SjjCjt \E\j\ 
S„C,r\E\.\ 


( 20 ) 


(19) 
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where 


j,t, !?= 1, 2 ,...i 


where 


S„C' jr (24) 

S„C'„~ \E*\ 

j= 1, 2, . . . i 

r=i+l, ...» 


Equations (17a) and (21a) state that the i-f-1 to n rows of 
the control matrix contain, only diagonal elements. Equa- 
tion (19) gives the required ratios between any two elements 
of any column of the O matrix. Equation (24) gives the 
required ratios of any element in any column of the C 
matrix to the diagonal element. 

Response equations. — The use of the operational functions 
to characterize each unit in the system allows algebraic 
representations of the relations among the variables of the 
system. With complete noninteraction specified, the system 
matrix of figure 5 becomes a diagonal matrix and each circuit 
acts as follows: 

y i =T,E Jk S kt C tj (Y J — Ljjyj) -f- S Ejjc V t (25) 

i=l t=l 

and 

x r =S rr C'„(X r -L„x r )+V r (26) 

When these equations are solved for the controlled variables 
as a function of the settings and the transient disturbances, 
the following equations are obtained: 


Vr 


S EjtSttCts 

*=i 


S (EjtSt&CkjLjj ) + 1 

*=i 


S E,*V k 


T,(E ilt S kk C ti L ii )+ 1 

t=i 


x r = 


S„C ' 


S TT C' TT L rT -\- 1 


■x t - 


Vr 


SrrC'rrLr, 


(27) 


(28) 


If the operational response ftmction of y } to Y s is defined as 
R if and the operational response function of x r to X T as R'„, 
equations (27) and (28) become 


n 


Vi—Rjr Yj—^t 1 ) Eji Vt 

t= i 


(29) 


where 


and 


x T =R'„-X-(R' rT L r r ~l)V r 


Bis 


S EjtSttCt] 

fc-i 

S ( Ejt Sue C k j Ljj) -j- 1 


R'r 


S„C' rr 

- S rr C' rr L rr +l 


(30) 

( 31 ) 


( 32 ) 


587 


The response of any y } to Y t insofar as the control matrix 
is concerned depends only on the elements in the jih. column 
of the control matrix (fig. 3). Inasmuch as for complete 
noninteraction only the ratios of the elements of that 
column of the control matrix are specified, the response can 
be independently set by choosing any one element of that 
column. If the conditions for noninteraction are used, 
Rj } is a function of any C. } 




(33) 


Equations (32) and (33) show the response of any controlled 
variable to its setting as determined by the control functions. 
The response of a controlled engine-independent variable 
does not depend on engine characteristics (equation (32)). 
The response of a controlled engine-dependent variable 
depends on the square engine matrix E* (fig. 2). 

An important result is obtained when equations similar 
to (29) and (30) are written for the errors in controlled 
quantities as follows: 

Y 3 -y s =(X-R ii L„) ( Y } -L s} g E*V^ (34) 


X r -X t =( 1 -R'rrEr) (X~ Xr V T ) (35) 


The errors in controlled quantities respond in the same 
manner to both setting and transient disturbances and the 
control configuration chosen can act both as a control to 
setting disturbances and a regulator to transient disturbances, 
at least with the transient disturbances assumed. 

Ideal control functions. — The control-design problem is one 
of fitting a control system to a particular plant or engine 
in order to control certain variables in some specified manner. 
The nature of the final system is characterized by the desired 
response of controlled quantities to outside disturbances. 
These responses may be exactly specified or given in general 
qualitative terms. The nature of the engine or plant to be 
controlled must be known in its static and dynamic states. 
Any method of explicitly solving for the control functions 
for any given plant and corresponding desired response 
action would give an exact answer to this problem. 

The use of the operational functions to describe the engine, 
control, and response functions and the algebraic methods 
employed herein lead directly to an explicit solution for the 
ideal control functions. Equations (32) and (33), which 
give the response actions in terms of the control function, 
can be solved for the control functions in terms of the 
response functions as follows: 


q n Bjj \E*j,\ 
" "“a-XXn) \E*\ 


(36) 


S„C„= 


R' 


1 — R' rr L n 


(37) 
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Choosing any desired response characteristics and solving 
for the required control functions to achieve this response 
are therefore theoretically possible. The use of these 
equations as well as the noninteraction conditions previously 
given for control design has certain limitations that pre- 
clude the use of this method in its most general sense. These 
limitations arise from the inability to incorporate exactly 
in a control the engine characteristics that the control is to 
match. A simple first-order response can be called for in 
equation (36) in spite of the higher-order transfer functions 
of engine and servo. This type of response is possible only 
if the control function can reduce the order of the transfer 
function by allowing factoring of the numerator and the 
denominator of equation (33) into exactly the same terms. 

The control functions found from equations (36) and (37) 
for responses that are stable would be expected to give 
satisfactory results for small changes in characteristics and 
in the engine characteristics the control function is to 
match. Special cases of engine characteristics may require 
that an unstabilizing term in the engine functions be 
factored out by the control. In this case, variations in 
engine or control characteristics would unstabilize the 
system. In such cases, equations (36) and (37) should be 
used with the proviso that the control does not decrease the 
order of the response function and the response functions 
must be chosen accordingly. 

The algebraic methods used in dealing with operational 
functions can be further generalized in terms of equivalent 
relations when the engine, the control, and the other unit 
characteristics are expressed in other forms. For instance, 
if the frequency-response characteristics are graphically 
known from tests or analysis, all the results described herein 
can be graphically performed in terms of the frequency- 
response characteristics in place of the operational functions. 

Summary of basic theory and results. — For an engine 
having n engine-independent variables and for which i 
engine-dependent variables are to be set and controlled (i ^ n) 
and for which n—i engine-independent variables are to be set 
and controlled, engine, control system, instrument character- 
istics, servo, and transient disturbances are 
Engine: 

yi=J±E ik -x k (8) 

*=I 

where 

j= 1> 2, . . .i 


Control: 


£*=SCit»-(i'»— y.)+ S G'tt-XXf — arj (10) 

V=\ /l“ t+1 


where 

Instruments: 


k= 1, 2 , ... 


3/ u — I'm ‘ Vv 

L/ift * Xp 


where 


®=1, 2, . . . i 


= ... a 


Servo and transient disturbance: 


z*= &**£*+ F* (12) 

where 

k= 1, 2, ... n 

The necessary and sufficient- conditions for aiiy'sctting J’ t 
to affect only y, are 

0 (17a) 

for all ji from f-j-1 to n, and 


SnC it \E* ts \ 
S„C H E*„\ 


(19) 


The necessary and sufficient conditions for any setting A”, 
to affect only x T are 


CV= o 

for all n from i+1 to n when 

sc , -T.\E*„\E tr 
S„C' r r \E*\ 


(21a) 


(24) 


For complete noninteraction, equations (17) and (19) must 
be true for all values of t from 1 to i and equations (21) and 
(24) must be true for all values of r from i+l to n. 

With complete noninteraction specified, the response of 
the controlled variables and their measured error signals to 
setting and transient disturbancesjirc 


and 


Yj-JZ l)E Jt Vt (29) 

X, = Ii'rr • A 1) Fr (30) 

Yi—yt=(\ —RtjLjj) ( F,-£„g E Jt F*) (34) 


where 


Xr-x r = (l -Ii\ r L„) (X r -L„ F r ) 


(35) 




^EjEftSktCtj 


}}— ~i 


\E*\S t ,C tj 


g(WA )+l 

(31) and (33) 


and 


R'> 


SrrC’ 


S„C'„L rr +l 


(ID 


(32) 
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With complete noninteraction specified, the explicit ex- 
pressions for the ideal control functions in terms of the de- 
sired response functions are 


op. Rn \E*j,\ 
" "“(1 -R«L it ) |E*| 

(36) 

O >0/ B rr 

” ” X _ R ’ ^ 

(37) 

SPECIFIC EXAMPLE 


If the turbine-propeller engine, which has the two inde- 
pendent variables fuel flow and blade angle, is considered, 
the defining differential equations obtained by linearizing the 
general functional forms of engine torque as a function of 
speed and temperature, propeller torque as a function of 
speed and blade angle, and temperature as a function of 
speed and fuel flow are 

A T + rA r = — bv)f 

(38) 

T=cii'j — eN 

(39) 


and from equation (24) for complete noninteraction, 

SnC'n ~\E* ll \E li E a b 

S 2i C' 32 \E*\ E n a 

where the determinant of a single-element matrix written as 
\E*\=\E* n \E u allows the use of the equations employed 
herein for this case. 

The control-system equations thus become 

P = C U (N, —N) +~ f? <r n <«/..-»/) 

— a On 

w r =C' 2i (w f ' l —w f ) 

The responses of speed and fuel flow to their settings 
become, from equation (33), 

x> f-E^j&iCii _ EuSuCu 

“ n -\E*\SnC n Ln + \E* n rE n S n CnLn+l 

aSuCn 

aS n C n L n — (1 rp ) 


When these equations are placed in operational form and 
speed and temperature are solved for in terms of fuel flow and 
blade angle, 


where 


N ~l + rp^ + l + rp- a ' 

(40) 

T- “ .0 + c(l-7+r3.) 

\-\-rp 1 + rp r 

(41) 



CASE i 


and from equation (32), 


22 SaP'nLu+l 

The expressions for the ideal control functions then become, 
from equation (36), 

o p Ru Ru 1 __ — Bn (l~r rp) 

1 — RnLii \E*\ ( 1 — En ( 1 — RnLii) & 


and from equation (37), 

SnC't 


1 


R' 22 

— R 22^2 2 


Controlling N and rc f . — For the case where n= 2 and i= 1, 
the engine matrix of figure 2 is obtained from equation (40) 
and 


En 1 + rp 

^=TfT 
1 +rp 

Xi — fS 

II 

s* 

< 

II 



The general control system is ' 

CniN, — iV)-J-C' / ia(w/,* — Wf) 

Wj— C u (N t —A 7 ) + C" 2 2(w/. , — E,) 

From equation (17a), the conditions for complete noninter- 
action give 

C 21 =0 


case n 

Controlling N and T . — For the case where n = 2 and i=2, 
the engine matrix is obtained from equations (40) and (41) 
and 


F - 

— a 


b 


~1 + rp 

12“ 

1 + rp 

F - 

at 


c(l~y+rp) 

■£*21- 

1 + rp 

■*-*22 

1 + rp 


X!=(3 x 2 =Wf 

yi~N y z =T 


The general control system is 

£= C n (N-N) + C U {T-T) 

m= C n {N-N)+ C 22 (T-T) 
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The conditions for complete— noninteraction give, from 
equation (19), 


(1— y+rp) 

b 


SnC n 


1 

fell 

1 

__6_ 

S22C21 

e\ 2 \~ 

E2I 

ay 


S 11 C 12 \E* 2 i 

L Zu 

S 22 C 22 j E*22 

! 


The control-system equations thus become 

^Cn(N-N)+~ f- 2 Cn(T.-T) 
— a On 


ay 




~ T ~ 6(1-7 +p)Sn 
The responses of speed and fuel flow to their settings become 

R n - 


\E*\8 n Cn 


SnCn 


a 


or 


-\E*\SnCnL n +\E*n\ (1-7+rj?) 

S 22 C 11 


R11- 


S 22 C 2l ill +T 


R’y 


\E*\S 22 C 22 _ SnCj 


28 


\E*\S i 2C22L 2 2~ j t-}E*22\ <y p j , 1 

022O22 Ar 2 aTT 


The expressions for the ideal control functions then be- 
come 

CY p Ru !- 2 '*nl —Rn (1 — 7 + rp) 

^^-l-BnLnWr^l-RnLn a 


S 22 C 22 - 


R‘. 


22 


E 


R* 


1 — R 22 L 22 I E* j (1 — R 22 L 22 ) c 


SUMMARY OF RESULTS 

The control of gas-turbine engines was limited to the 
linear assumptions for the engine, the control system, and 
the other units in the system. Operational functions were 
used for these units, the system was generalized to allow 
for complex engine types having any degree of freedom, and 
algebraic methods were employed in setting up the char- 
acteristics of the engine, the control, the servo, and the 
measuring instruments. These algebraic methods were em- 
ployed to combine the various separate units into a system 
showing how various effects are transmitted throughout 
the entire system. A condensed matrix form was used 
to illustrate the principal effects on controlled variables. 

By use of either the index or matrix form of the algebraic 
equations of the system, the following principal conditions 
were explicitly expressed: 

(a) The conditions that must be imposed on the 
control functions so that any setting disturbance will 
affect only the variable being set 

(b) The conditions that must be imposed on the 
control functions for complete noninteraction so that 
every setting disturbance will affect only its corre- 
sponding variable 

(c) The response functions, which characterize the 
controlled-system response action, in terms of the 
control-system functions 

(d) The ideal control-system characteristics in terms 
of any desired response action. 

A specific example using the turbine-propeller engine 
having two degrees of freedom was given to illustrate the use 
of the previous results. 

Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, April 25, 1949. 



APPENDIX 

CALCULATION CONDITIONS 


Necessary condition. — The general equations for the re- 
sponse of the controlled variables are given by equations 
(Al) and (A2) (equations (14) and (15), respectively, in the 
text) as follows: 


CONDITIONS ON C MATRIX 

where 


Vi — Xj EjtSatCt 0 ( F, — L„y r ) -|- 

t=l v= 1 


S S Ejk&up' t/JCf — + 2 E jk Ffc 


(Al) 


where 


3—1, 2, . . . i 


x*^'£,S*C k ,(Y,-L„y,)+ S 

« = 1 n^t-rl 


where 


(A2) 


k= 1, 2, ... n 


If the effects of any setting Y t are separated, equations 
(Al) and (A2) become 


Y,—L, a y,) -f- 

S S EjtSuC'tpiXp — LppX^) -f- 
*=I ,.=£+! 

’^jEj k V k -\-'^^E jk SichCict(Yi — L u y t ) (A3) 

x k =’^ J S kk C^(J , — L„y,)-\- S/^C'tpiXp — LppX ^) + 

tt=* I 

I a*{~ S*Cu{ 1 1 — L tt yt) (A4) 

If any setting F t is to affect only y t , then any y s is not affected 
by Y t (jy^t) and any controlled x t (&=i+l to n) is not af- 
fected by Y t . If equation (A3) is considered for any j t, 
then the only term that can be a function of l 7 ! is the last 
one. Because of the arbitrary nature of the instrument 
function, Y t ^L tt yt so that 


where 


^EjiSuC&t — 0 




(A5) 


If equation (A4) is considered for k=i-\-l to n, the only 
term that can be a function of F, is the last one; therefore 


SaCtt—O 

Cju — o 


(A6) 

(A7) 


Jr=i+1, ... 7i 

If equation (A7) is used, equation (A5) becomes 


^jExS^Cu — 'E'djtEACtt 
k -1 k=l 


where 

and 

when 

and 

when 


3 — 1 , 2, . . . i 

8]t=0 

S Jt = 1 
j=t 


The solution of this set of i equations in the i unknown 
S kk O kt gives 

i,±,\E\ i \8, t E lk S kk C tt 

Cl /~1 P=I k -1 




and because 

ll 

o 

when 

P7*t 


S pt =l 

when 

p=t 


S]]C jt = 


\E*u[ 2 EaS^Ckt 

i 


where 


i=l,2, . . . i 


(AS) 


Because j is any number (1,2, ... i), a similar expression 
for any other S„ G, t is 


S„C. t =- 


\E* t t \ S EttS&Ctt 


k=l 


The ratio of equations (A8) and (A9) becomes 


where 


s»c Jt 

\E* ti 

s„c, t 

\E* t . 

II 

h- 4 

JO 

. . . i 

<5 

II 

h-* 

JO 

. . . i 


(A9) 


(A10) 


Equations (A7) and (AlO) thus follow necessarily from the 
noninteraction condition specified. 
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Sufficient condition. — If equation (A10) holds, 


where 


/ o ft 

•STA T7I £>!}'-' it 1=1 

Z-J E'TS a p I TT* r 

j=l O n O vt !■& «1>I 


r= 1,2, . . . i 

From the theory of determinants, the expression 


when 

when 

Therefore, 

when 


-£E fi \E* tJ \=0 

j=i 

r^t 

±E rj \E* ti \=\E*\ 

= 1 

r=t 


y 2 E ri s„c i( =(j- 

j = i 


(All) 


TjZt 


If equation (A7) also holds, then equation (All) becomes 
±E ri S fi C it =0 (A 12) 


where 


r?*t 

#•=1,2, . . . i 


where 


i=l,2, . . . i 


If equation (A14) is considered for any £(i+l to n) where 
k^r, the only term that can be a function of A’ r is the last 
one and therefore 

Su.C’^0 (A 10) 

and 

<7V=* 0 (A 17) 

where 

k 9 * r 

If equation (A17) is used, equation (A15) becomes 


where 


2 EfrSutC'kr = — Ejr S„C'r 

1=1 


3- 1, 2, . . . i 


The solution of this set of i equations in the i unknown 
SkkC’tr gives 

-’h\E\ f \E pr S„C ) rr 


S ru p-i 

U'-' Jr 


or 


m 

of fit SI E*„ } \E pr 

ir _ l 


S„C' 


where 


W\ 

j— 1. 2 .i 


(A 18) 


If equations (A12) and (A7) hold, then any setting Y, will 
affect only y t (fig. 5). Equations (A7) and (AlO) thus are 
sufficient for the noninteraction condition specified. 

CONDITIONS ON C' MATRIX 

Necessary condition. — When the general equations for the 
response of the controlled variables (equations (Al) and 
(A 2)) are written separating the effects of any setting X„ 
these equations become 

Vj— 2) S E Jk SkkCto(. Y, — LcvVt) 4" 

1 i>=l 

2 2 E jk S kk C' iiiXp — L^Zp ) + 

k= 1 p=i-H 
n*r 

p^Vt+gEAC'UX-LrrX,) (A 13) 

%k—'^jSi ! kCi C o(Y c —L sv y,)-f 2 SmC'i-JX^ — L^ z ^) + 

f=l #i=f+l 

Vt+StoC'UXr-LrrX,) (A 14) 

If any setting X r is to affect only x T , then any y } is not 
affected by X r and any controlled x k (k = i+l to n) is not 
affected by X r for k 9 ^ r. In equation (A13) the only term 
that can be a function of X r is the last one. Because of the 
arbitrary nature of the instrument function, X T 9 * L„x, and 

S E^ SttC'tr =0 (A 1 5) 

k= 1 


Equations (A17) and (A18) thus follow necessarily from the 
noninteraction condition specified. 

Sufficient condition. — If equation (A18) holds, 

f rp n ft/ zL )S E v f\E pf\E nr 

^sT\ J-'vJ&JJ'-' Jr ■ p=lj-l « 

U SrrC’rr + \E*\ ° 

From tlilTtheory of determinants applied to the second term 

/-I O rr O rr 

and 

2 E,fS„C'jr+E„S rr C t r =0 

j=i 

If equation (A 17) also holds, tlie previous equation becomes 

2 E'fSffC'ir^O (A 1 9) 

1 

If equations (A17) and (A19) hold, then any setting X, 
will affect only x r (fig. 5). Equations (A 17) and (A 1 8) thus 
are sufficient for the specified noninteraction conditions 
to hold. 
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